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Abstract 


In  this  report  entitled,  "Refined  Test  Methods  In  Evaluating  Brittle 
Material  Fracture  Strengths,"  Contract  No.  N00019-79-C-057i» ,  is  a  follow¬ 
up  of  the  work  done  under  NAVAIR  No.  N00019-78-C-0520,  entitled,  "Fracture 
Prediction  in  Brittle  Materials.". 

In  the  former  contract,  an  analytical  system  was  devised  to  generate 
fracture  data  based  on:  (1)  conversion  of  triaxial  fracture  data  to  a 
single  variable,  $,  the  angle  on  the  fracture  surface  from  point  of 

A 

fracture  to  the  n*(l,l,l)  direction,  (2)  the  creation  of  a  simple 
test  procedure  and  specimen  design  to  find  a  variety  of  points  on  the 
fracture  surface,  and  (3)  the  design  of  experiments  to  employ  this  testing 
system.  In  the  first  contract,  both  alumina  and  mulllte  specimens  were 
tested  in  four-point  loading  and  these  results  were  reduced  to  a  Weibull 
distribution  by  determination  of  the  two-parameter  coefficients. 

In  this  contract,  a  torsion  test  method  and  apparatus  were  designed 
and  implemented.  Both  alumina  and  mullite  specimens,  the  same  geometry 
as  was  used  in  the  four-point  bend  test  of  the  previous  contract 
were  loaded  in  torsion  to  failure  at  another  angle  on  the  fracture 
surface,  or  at  $-90  to  the  n* (1,1,1)  direction.  The  results  are  given 
for  alumina  and  mullite  tests. 

Finally,  a  general  test  utilizing  'thick'  disks  of  each  material 
is  described.  The  literature  seems  to  lack  a  simple  reliable  disk  test 
for  which  a  simple  stress  distribution  is  readily  available.  The 
test  produces  an  axially  symmetric  state  of  stress  but  one  that  embodies 
a  variety  of  orientations  to  the  n- (1,1,1)  direction.  The  use  of  a 
hydraulic  ram  gives  rise  to  an  easily  predictable  stress  field.  The 


expected  Weibullian  scatter  was  detected  in  both  alumina  and  mullite 
specimens  which  were  tested  in  both  as-cast  and  ground  conditions. 

No  problems  attributable  to  lack  of  flatness  in  these  low-strain-to- 
f allure  materials  were  detected  In  as-cast  samples.  The  alumina 
samples  tended  not  to  show  fracture  propagation  throughout  the  thick¬ 
ness  while  the  mullite  samples  did  demonstrate  this  behavior.  Because 
it  was  sometimes  difficult  to  determine  precisely  what  pressure 
produced  incipient  cracking  in  alumina,  a  differential  current  detector 
employing  a  conductive  brittle  coating  has  been  proposed — the  coating 
brittleness  being  chosen  to  match  the  specimen  strain-to-fallure  . 
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Chapter  1 


THEORY  OF  THE  TORSION  TEST 

A  circular  prism,  rod.  or  tube  can  be  subjected  to  a  state  of  pure 
torsion  along  its  gage  length.  Fig.  1.1. 


Figure  1.1  Biaxial  Stress  Envelope  for  Circular  Cross-Section 
Tube  Subjected  to  Pure  Torsion 

The  Biaxial  stress  state  depicted  by  Mohr's  circle.  Fig.  1.2.  is 
applicable  to  every  point  of  the  specimen  if  edge  effects  are  neglected. 
As  demonstrated  by  the  results  in  Chapter  3.  the  edge  effects  did  not 
contribute  to  testing  error  in  a  significant  way.  However,  almost  some 
stress  concentration  is  to  be  expected  in  any  experiment  and  attempts 
should  be  devised  to  minimize  it  as  with  the  special  end  caps  and  ad¬ 
hesives  described  in  Chapter  2. 


Figure  1.2  Mohr's  Circle  in  Pure  Torsion 
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The  probability  of  failure  is: 


F(T)  -  1  -  exp (-Bn) 


1 .1 


The  risk  of  fracture  or  rupture  can  be  expressed  easily  In  integral 
form  since  the  radial  volume  element  and  stress  level  are  simply  related. 
Further,  a  scalar  variable,  T,  that  is  one  possible  contraction  of  the 
stress  state  torsion  is  used  in  the  integral  to  obtain: 


Br.  =  Lc 


//  (¥■  ■  °“)m  rd9 


dr 


where  c  = 


Vun«V 


m  =  Tr.  and  L  -  the  gage  length. 


1 .2 


End  Caps  filled  with  Adhesive 


Ri 


Fig.  1.3  Pure  Torsion  for  a  Circular  Cross-Section  Tube 


Using  the  same  approach  as  in  IU  rods  and  tubes  are  separately 
considered  : 


Two  cases  occur; 


(a)  0^?.  ,  C  Ro  ; 

tn 


c  J 
Uu 

where  R^.  “  -  ,  and  B  becomes 

th  x  n 


21TLC 


Br.  = 


»  rr*  \ 

lD+1,(j) 


/T  r  \wl 

r(r  •°“j 


(¥  •  -)' 
(!) 


ir+2  i 


Cm+2) 


Ro 


Rth 


(b)  Ro^?.^  ;  In  this  case; 

Bn  «  0 . 


(2)  Tube  (Ri  ^  0) :  Three  cases  occur; 

(a)  Ro>Rth>Ri 

Bn  Is  expressed  by  Eq.  (1.3) 

(b)  0_<Rth_<Ri 

Bn  is  expressed  by  Eq.  (1.3)  except  the  lower  limit 
of  integration  is  replaced  by  Ri. 

(c)  Rth>Ro;  in  this  case; 

Bn  -  0., 

Equation  (1.3)  has  the  same  form  as  the  case  of  the  uniaxial  tension 
field  and  can  be  correlated  with  either  the  three-parameter  Weibull  distri¬ 
bution  or  the  two-parameter  distribution  in  an  iterative  programming  scheme 

similar  to  UNIAX. FOR.  [1] 

Any  of  the  scalars,  the  shear  stress,  T;  the  effective  tensile 
principal  stress,  Oj«  T  r/J;  or  the'biaxial  intensity o,  can  be  used  in 

Eq.  d.3)  for  the  study.  This  intensity  might  be  developed 
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in  several  ways  including  the  correlation  of  test  results  with 
1)  the  most  tensile  component  of  the  stress  tensor  or  2)  the  strain 
energy  of  distortion  [2]. 

In  either  case,  each  element  is  considered  to  be  under  a  torsional 
stress,  where  t  =  Tr/J  is  assumed,  and  the  risk  of  rupture  for  the  whole 
specimen  is  obtained  by  integrating  the  risk  of  rupture  of  each  infinites- 
inal  element  over  the  gage  volume  of  the  specimen.  Thus,  1.3  applies. 
Under  this  pure  torsion,  the  stress  field  xrz=Tr/J  results  from 

which  the  principal  direction  stresses  are: 

o,  =  t  ;  c  «-t  1 .5 

1  rz  2  rz 

Normally  the  compressive  component,  a would  be  neglected  and  a 
uniaxial  stress  state,  ,  would  be  assumed.  In  this  case  o=  o^= Tr/J. 

If  distortion  energy  is  employed  to  define  then  under  pure  torsion: 

u’ .  (w2  ♦  (v°3>2  *  ‘W2 


For  this  study, compression  was  ignored  and  o,  =  t  =  Tr/J 


led  to  1.3- 
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Chapter  2 

DESCRIPTION  OF  THE  TORSION  TEST 

The  torsion  test  method  required  the  creation  of  a  torsion  test 
loading  system  designed  to  operate  In  conjunction  with  a  standing 
universal  testing  machine,  Figure  2-1.  The  actual  torsion  test  fixture 
designed  to  accept  a  specimen  between  two  universal  joints.  Figure  2.2, 
was  also  designed  to  eliminate  all  loads  except  pure  torsion. 

The  pictorial  view.  Figure  2*3,  demonstrates  the  current  configu¬ 
ration  in  general  while  the  component  dimensioned  drawings  provided 
detail.  Figure  2-4,  5.  The  experiment  was  first  conceived  as  one  that  did 
not  require  a  locking  bar  behind  each  collet.  It  was  found,  however, 
that  the  experiment  was  unsuccessful  because  the  hardened  collets  or 
chucks  slipped  before  fracture  stress  was  obtained.  Naturally,  these 
early  tests  were  conducted  on  alumina,  the  stronger  of  the  two  materials. 
The  original  sleeves.  Figure  1*3,  were  thus  machined  to  allow  for  a  slot 
and  the  experiment  repeated. 

The  next  problem  was  one  of  adhesives.  A  rough  calculation  indi¬ 
cated  that  an  elastomer  adhesive  would  need  to  withstand  up  to  3000  p.s.l. 
in  shear  to  sustain  a  load  sufficiently  high  to  fracture  alumina.  After 
consultations  with  3M  adhesive  experts  a  film  adhesive  —  3M 
AF-126-2  —  was  selected  .  Surfaces  of  both  the  ceramic  and  metal  sleeve 
were  prepared  carefully  to  permit  good  bonding.  The  adhesive  was  cured  as 
specified  and  several  specimens  were  tested. 

The  results  of  these  tests  were  disappointing  and  all  joints  failed 
at  no  greater  than  25Z  of  the  expected  alumina  fracture  torque.  Perhaps 
this  problem  was  caused  by  (1)  use'  of  a  ’bad'  or  degraded  batch  of  adhesive 
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Figure  2.4  Torsion  Testing  Fixture-Di 
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Figure  2,5  Detail  of  Toralon  Spa clean  End  Grlpa 


as  affected  by  transportation,  or  (2)  by  the  lack  of  sufficient  interference 
pressure  during  the  curing  process.  The  only  available  pressure  during 
curing  was  obtained  through  the  interfemece  fit  between  the  sleeve  and 
alumina  specimens  largely  due  to  thermal  expansion  of  the  adhesive. 

The  second  adhesive  tried  was  3M  Structural  Epoxy,  Scotch-Weld  Brand, 
No.  2216  B/A.  This  gray  two-part  epoxy  had  less  advertised  strength,  but 
it  proved  far  more  suitable  even  though  a  room  temperature  cure  was  used. 
Only  in  a  very  few  cases  did  this  adhesive  fail.  In  all  cases  the  epoxy 
was  burned  out  after  each  test  and  a  new  specimen  was  bonded  again.  The 
burnout  of  end  caps  meant  that  fewer  caps  needed  to  be  machined,  thus, 
reducing  expense  considerably.  The  epoxy  was  also  used  to  ‘pot*  the 
ceramic  disks  in  place  on  copper  gaskets  in  the  disk  tests  as  described 
in  chapters  A  and  5. 

The  ceramic  specimens  used  were  made  by  McDanel  Refractory  of 
Beaver  Falls.  Pennsylvania.  Their  high  density  988-Alumina  and 
MV33-Mullite  ceramic  tubes  were  the  specimens  for  both  the  bending 
tests  of  [1]  and  these  torsion  tests.  The  same  materials  and 
batches  were  also  used  for  the  disk  tests  of  this  report.  For  the 
disks,  the  method  of  casting  was  of  necessity  different.  The  tubes 
were  extruded.  Table  2.1  describes  tube  geometries. 


Alumina  (998) 

Mullite  (MV33) 

1  “ 

Mean  and  SD  of  Maximum 

OD 

1 .0010+  .0021 

1 .0290+  .0065 

Mean  and  SD  of  Director 

1  to  Max.  OD  Direction 

.9985+  .0027 

1  .0236+  .0082 

Mean  and  SD  of  Minimum 

ID 

.6156+  .0203 

.7276+  .0087 

Mean  and  SD  of  Direction 

1  to  Min.  ID  Direction 

.6209+  .0234 

.7367+  .0074 

Estimated  Outer  Radius 

.4999 

.5132 

Estimated  Inner  Radius 

.3091 

.3661 

Test  Length 

BH 

7.3 

Overall  Length 

11  .0 

11  .0 

Number  of  Samples  Used  for 
These  Statistics 


50 


50 


Chapter  3 

RESULTS  OF  THE  TORSION  TESTS 

Alumina  998  Results 

Results  for  the  torsion  test  of  alumina  are  listed  below  in  Table  3-1. 

It  should  be  noted  that  most  entries  were  considered  successfully  tested. 
Entries  ranked  21  and  22  were  used  as  'good  tests'  even  though  they  did  not 
fracture;  their  epoxy  bond  failed  first,  but  had  a  very  high  stress  level. 

For  that  reason  and  because  the  high  end  of  the  curve  might  be  more  accurately 
depicted,  it  was  decided  to  include  them.  Each  of  these  other  specimens 
actually  fractured  in  torsion  and  their  end  caps  were  removed  by  heating 
the  failed  specimen  end  cap  with  an  acetylene  torch  to  burn  out  the  adhesive. 
The  cap  was  cleaned  in  a  solvent,  sandpapered,  then  reused.  Thus,  six 
sets  of  caps  were  required  so  that  three  were  bonded  to  specimens  being 
tested  while  three  more  were  being  cleaned  on  any  work  period.  At  two 
points,  the  mild,  cold-rolled  steel  caps  became  deformed  in  the  slot  through 
Dlastic  deformation  so  that  new  caps  had  to  be  manufactured,  c.f..  Fig.  2.5. 

Table  3.2  lists  all  tests  on  alumina  that  were  rejected.  In  the 
great  majority  of  such  tests  the  specimens  slipped  in  the  grips  indicating 
bond  failure.  Two  specimens  that  originally  slipped  at  low  levels  were 
retested  with  new  adhesive  and  were  successfully  tested  to  failure.  Most 
tests  in  which  bond  failure  occurred  were  never  retested  successfully. 

This  is  thought  to  be  because  those  torsion  specimens  were  found  to  be 
slightly  undersized .  thus,  contributing  to  a  thick  bond  of  reduced  sheer 
strength . 

It  was  decided  that  any  test  which  did  not  demonstrate  failure  through 
the  specimen  center  would  be  rejected.  This  is  because  tne  most  accurate 
representation  of  the  Weibull  Parameters  was  desired.  Flaws  near  the 
attachments  might  Dropagate  to  failure  by  bond-induced  stress  concentrations, 
thus  invalidation  the  results.  No  such  examcles  occured:  however,  it 
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should  be  pointed  out  that  most  specimens  exhibited  a  spiral  shaped 
crack  through  the  grip  attached  to  the  torsion  arm  on  the  testing 
machine  -  the  active  grip.  This  phenomenon  was  thought  to  result  from 
the  release  of  strain  energy  in  the  heavy  assembly  through  springback. 

All  such  specimens  failed  in  the  gage  -  length  end  active  grip. 

The  data  were  analyzed  as  a  two-parameter  Weibull  distribution 
through  the  use  of  an  iterative  programming  code  in  Appendix  I. 

AUTO. FOR  and  AGARWA. FOR.  The  results  for  the  alumina  samples  are 
listed  below: 

Two-Parameter  Family  for  Alumina  998 

0.617E-28  (in.2m  lbs."m):  2.708E-52  (M2m  N“m) 

6.085  ''dim  less):  6.085  (dim  less) 

0  'psi):  0  'MPa) 

4.325  E+4  (psi):  292.2  (MPa) 

Residual  error  from  curve  fitting,  Re  =  .  33^0E-O1 

N  *\.  2 

where  Re  =  l  W.  «  P*) 

1=1 

N  =  Number  of  entries 

p  =  Probability  of  Fracture  used  as  input 

?  =  Derived  probability  of  Fracture.  Table  3-3 

=  Weight  of  each  Data  Entry. 

Because  probabilities  of  fracture  vs.  load.  Table  3.1>  were  generated 
by  the  rank  method,  the  use  of  the  three  parameter  distribution  was  not  used. 
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Taken  as  same  rank  with  twice  normal  weight 

These  tests  were  included  even  though  there  wap  bond  failure  rather  than 
fracture  because  the  load  at  failure  was  considerably  above  the  mean  of  other  tests 


Table  3-2.  Data  On  All  Other  Alumina  Tests 


Specimen 

Number 

Torgue  at  End  of  Test 
(Foot  Founds) 

Remarks 

ATI 

132.50 

No  end  cap  used  —  initial  set-up 
test. 

AT2 

100 

Slipped  in  grip  of  initial 
design  end  slot. 

AT  3 

100 

Slipped  in  grip  of  initial 
design  without  end  slot. 

AT4 

259.38 

First  modification  of  Initial 
grip  with  no  end  slot  but  with 
roughened  surface  in  contact 
with  collet;  it  slipped. 

AT5 

259.38 

First  modification  of  initial  grip 
with  no  end  slot  but  with 
roughened  surface  in  contact  with 
collet;  it  slipped. 

AT22 

406.09 

The  following  tests  slipped  in 
the  bond  with  grip  of  latest 

AT53 

200 

design  —  some  were  retested 
and  still  slipped.  The 

AT54 

200 

specimens  were  discarded  because 
they  exhibited  smaller  than 

AT55 

200 

average  outer  diameters  requir¬ 
ing  a  thick  bond  line;  thus. 

AT63 

200 

making  for  a  weakened  test  with¬ 
out  machining  undersized  end 
caps. 

The  actual  predicted  fracture  probability  vs.  load  and  input 
fracture  probability  vs.  load,  are  given  in  Table  3.3  below.  The 
probability  of  failure  for  a  1  inch^  specimen  is  given  by  Graph 
3.1  for  equation  parameters  3-1. 


Specimen 

Rank 

Probability  of 
Fracture  used  as 
Input 

Probability  of 
Fracture  from 
Equation  3.1 

Weight 

1 

.04348 

.01997 

1 

2 

.08696 

.06428 

1 

3 

.13044 

.13205 

1 

4 

.17391 

.15447 

1 

5 

.21739 

.16639 

1 

6 

.26090 

.22631 

1 

7 

.30435 

.27519 

1 

e 

.34783 

.40665 

1 

9 

.39130 

.43272 

2 

10 

.43478 

.45462 

1 

11 

.47826 

.52563 

1 

12 

.52174 

.53850 

1 

13 

.56522 

.56273 

1 

14 

.60870 

.62682 

1 

15 

.65217 

.64900 

1 

16 

.69565 

.66534 

1 

17 

.73913 

.70284 

1 

18 

.78261 

.72108 

1 

19 

.82609 

.78145 

1 

20 

.86956 

.89587 

1 

21 

.91304 

.99549 

1 

22 

.95652 

.99959 

1 

TABLE  3. 

3  Assumed  and  Determined  Probabilities  of 
Alumina  998  Original  Test  Data 

Fracture  : 
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A  photograph,  Figure  3.1,  demonstrates  the  fixture  system  immediately 
following  one  of  the  tests.  Note  the  lack  of  crack  propogation  into  the 
end-caps  and  the  spiral  nature  of  the  fjracture  surface'.  Further  examples 
of  Alumina  specimens  just  after  failure  depict  samples  AT60  and  AT62. 
Attention  again  is  drawn  to  the  lack  of  any  end  cap  fracture  zone  that  was 
detectable.  Again,  a  spiral  fracture  surface  is  evident  that  seems  to 
initiate  and  terminate  on  a  material  bifurcation  point  or  cusp  persumably 
representing  weakend  material  -  Figures  3.2  and  3.3.  Generally  the 
Alumina  tubes  exhibited  very  high  strengths  in  comparison  to  mullite. 

Some  specimens  did  not  fracture  until  surface  tractions  exceeded  50  ksi. 


Mullite  MV33  Results 

Mullite  Specimens  selected  randomly  for  the  torsion 

test  were  labeled  as  depicted  in  Table  3.^. The  Mullite 

(MV33 )  specimens  exhibited  somewhat  less  than  J  the  mean  strength  of  the 

Alumina  tubular  specimens  and  were  of  a  different  mean  dimension,  Table 
2.1.  Because  the  Mullite  was  so  much  weaker,  no  end  caps  or  adhesive 


bonds  failed  in  these  tests.  For  these  specimens,  the  Weibull  parameters 
were  found  to  be: 


c  =  .958E-27  (in2m  lbs"m) 

m  =  6.1  78  (dim'less ) 


1.848E  -51  'M2m  N"m) 


6.178  (dim'less) 
0  (MPa) 

o«  2 . 3628E+4  (o=  c_1 /m) (psi ) :  1629  (MPa) 


a=0  ( psi ) 
u 


o  o 

Residual  Error  from  curve  fitting  *  .3846E-01 

As  a  measure  of  the  fit  of  the  derived  distribution  with  the 

initial  data,  Table  3.5  provides  a  comparision.  The  graph  3-2  depicts 

the  shape  of  the  fracture  curve  for  a  1  inch"*  volume  subjected  to  uniaxial 

tension. 
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Figure  3.3.  Felled  Alumina  Torsion  Specimen  AT62 
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Table  3.4  .  Torsion  Experimental  Results  —  Mullite  MV33 
(Data  from  24  good  tests) 
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Probability  of 

Probability  of 

Fracture  used 

Fracture  from 

Specimen  Rank 

as  input 

Equation  3.2 

Weight 

1 

.04546 

.01949 

1 

2 

.09091 

.04616 

1 

3 

.13636 

.07059 

1 

4 

.18182 

.12260 

1 

5 

.22727 

.29362 

1 

6 

.27273 

.30784 

2 

7 

.31818 

.35139 

1 

8 

.36364 

.36270 

1 

9 

.40909 

.42203 

2 

10 

.45454 

.40869 

1 

11 

.50000 

.49429 

2 

12 

.54546 

.56249 

1 

13 

.59091 

.59187 

1 

14 

.63636 

.62112 

1 

15 

.68182 

.62407 

1 

16 

.72727 

.67853 

1 

17 

.77273 

.68475 

1 

18 

.81818 

.86166 

1 

19 

.86364 

.91127 

1 

20 

.90909 

.93219 

1 

21 

.95455 

.96981 

1 

Table  3.5. 

Assumed  and  Determined 

Probabilities  of 

Fracture 

Mullite  MV33  Original 

test  Data 

The  computer  programs  used  to  derive  these  data  are  interative 
in  nature  which  may  seem  inconsistent  with  the  existence  of  a 
closed  form  for  the  Weibull  volume  integral  (1.3)  and  simple 
Least  Square  solution  available  [1  for  circular  torsion  specimens. 
This  approach,  however,  allows  other  minimization  parameters  to  be 
used  as  well  as  a  3-parameter  model  as  appropriate. 


It  is  interesting  that  the  fracture  of  the  mullite  specimens  differ 
substantially  from  the  fracture  of  the  alumina  specinmns.  The  difference. 


one  of  number  and  size  of  the  fracture  fragments,  is  similar  to  that 

observed  in  the  bend  tests  although  on  a  more  dramatic  scale.  Figures  3-4, 
3.5,  3.6,  and  3.7  show  that  the  fracture  of  mullite  produces  a  great  many  more 
fragments  than  the  alumina;  the  uniformity  of  the  fracture  patterns  with 
generally  a  spiral  geometry  is  unmistakable.  This  uniformity  is  thought  to 
be  a  result  of  (1)  the  nearly  constant  strain  energy  density  achieved  in 
this  type  of  test,  and  (2)  the  apparent  lack  of  an  appreciable  stress  con¬ 
centration  at  the  caps.  Some  of  the  specimens  in  the  figures  have  failure 
well  into  the  caps.  Figure  3-4»  through  3.7,  and  some  do  not,  Figure  3.8. 


Figure  3-4.  Mullite  Failure  in  Torsion..  Specimen  MT87 


V 


Figure  3.5.  Mullite  Failure  in  Torsion,  Specimen  MT88 


X 


Figure  3.6.  Mullite  Failure  in  Torsion,  Specimen  MT89 
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Figure  3.1.  Mullite  Failure  in  Torsion.  Specimen  MT90 


The  existence  of  specific  bifurcation  crack  initiation  points  again 
<•>  is  in  some  evidence  in  the  Mullite  specimens.  This  material  demonstrates 

greater  brittleness  than  Alumina,  a  finding  also  to  be  noted  in  the 
pressure  tests  and  in  the  4-point  loading  tests  of  [1]. 


38 


Chapter  4 


Discription  of  the  Hydraulic  Test 

The  experimental  designs  and  methods  employed  for  the  tube  specimen 
tests  are  characterized  by  nearly  uniform  stress  states  in  the  test  section 
except  for  the  bending  test  where  a  constant  principal  stress  ratio  was 
established.  If  there  is  symmetry  of  Weibull  parameters  as  a  function  of  the 
direction,  4  the  angle  between  that  point  on  the  fracture  surface  and  the 
n  =  (1,1,1)  1  direction,  then  the  identification  of  these  parameters  at  a 

vr 

sufficient  number  of  values  of  4  defines  fracture.  The  table  below  describes 
values  of  4  obtained  by  common  tests . 


4  Degrees 

Test  Type 

54.7 

Uniaxial  tension  or  pure  bending 

39.2 

Internal  Pressure 

90.0 

Torsion 

125.3 

Compression 

Table  4.1.  Fracture  Surface  Angles  for  Typical  Types  of  Tests 

While  neither  compression  nor  internal  pressure  tests  were  conducted, 
they  could  be  employed  on  the  same  tube  specimens.  The  test  developed  and 
described  in  this  chapter  supplies  a  variety  of  stress  ratios  by  hydraulically 
loading  disk  specimens  simply  supported  along  their  edges,  Figure  4-1.  This 
test  is  necessary  to  complement  fracture  theory  data  gleaned  from  tube  specimen 
employing  a  single  4»  Independence  of  the  action  of  each  principal  direction 
stress,  a  commonly  applied  assumption,  can  be  evaluated  m  this  way  [4],  In 
addition  this  hydraulic  test  may  be  an  alternate  to  disk  bend  tests  that 
roughly  approximate  this  method  [5,6]. 


Figure  4.1.  Alumina  and  Mullite  Disk  Specimens  and  Loading  Conditions 


This  also  provides  an  additional  geometry  for  application  of  .'leibull  theory. 
Both  'as-cast'  samples  and  disks  machined  to  ensure  flat  parallel  surfaces  were 
used  as  specimens  in  both  Alumina  and  Mullite  from  the  same  batches  as  those  of 
the  tube  specimens.  D,  d  refer  to  the  radius  and  thickness  of  'as-cast'  disks 
and  D,  d  refer  to  those  measurements  of  machined  disks.  The  machined  specimens 
were  ground  for  flatness  and  for  the  removal  of  sharp  edges,  Fig  4-1.  resulting 
in  a  radius  of  about  .125".  There  was  concern  that  'as-cast'  disks  would  be  pre¬ 
disposed  to  warpage  thereby  causing  considerable  parasitic  stresses  during  test¬ 
ing.  In  order  to  examine  this  potential  problem,  'as-cast'  and  ground  specimens 
in  both  alumina  and  mullite  were  all  tested  to  failure.  Their  dimensions  and 
deviations  are  given  in  Table  4-2  in  reference  to  the  geometries  of  Figure  4-1. 


D 

d 

d 

Sample  Type 

Diameter  (in) 

Thickness  at 
Center  (in) 

Thickness  .25" 

From  Edge  (in) 

Alumina  -  'as-cast' 

2.02101.0140 

.25221.0038 

.2519*. 0034 

Alumina  -  ground 

2.026 3± . 016 1 

.27301.0060 

.27 36 ±.006 3 

Mullite  -  'as-cast' 

2.07211.0172 

.2507*. 0012 

.25081.0014 

Mullite  -  ground 

2. 0533*. 0196 

.24881.0058 

.25021.0048 

Table  4.2  Disk  Specimen  Dimensions  -  means  and  standard  deviations 


In  reference  to  Table  4.2,  it  apperars  that  the  supplier  selected  thicker  alumina 

samples  for  grinding  and  that  grinding  did  not  improve  the  overall  thichness  devi¬ 
ations. 


While  there  have  been  numerous  proposed  disk  specimens  for  brittle 
materials,  the  majority  of  such  tests  in  checking  Weibull  parameters  are 
of  little  value  because  either  (1)  they  produce  inordinately  complex 
stress  fields,  (2)  the  mode  of  testing  imparts  high  stress  over  a 
small  volume,  or  (3)  the  small  strain  to  failure  requires  the  use  of 
expensive  high-precision  samples  and  fixtures. 

The  test  fixture  of  this  report  consists  of  a  small  hand-operated 
hydraulic  pump  capable  of  10  ksi  pressures  and  a  pressure  gage  and  tubing 
attached  to  the  specimen  testing  assembly,  Figure  4.2.  In  order  to  provide 
for  simple  support  around  the  side  of  the  specimen  and  to  discourage  leak¬ 
age  of  hydraulic  fluid,  a  soft  copper  gasket,  cut  to  fit  the  assembly  bolt 
pattern  is  bonded  with  3M  structural  epoxy  to  the  edge  of  the  specimen. 
Figure  4  .3. 

The  specimen  with  copper  gasket  is  sandwiched  between  heavy  steel 
members  and  held  together  by  stiff  preloaded  bolts.  One  side  of  the 
specimen  is  in  contact  with  pressurized  hydraulic  fluid  while  the  other 
side  is  vented  to  the  atmosphere  through  a  1-1/2  inch  hole.  Figure  4-4. 

The  difference  in  a  diameter  of  the  specimen  and  vent  provides  a  lip  of 
1/4  inch  allowing  contact  with  the  soft  copper  seal  to  achieve  a  simply- 
supported  condition.  More  detail  is  provided  by  Figures  4-5  through  4-7. 

With  each  test  the  apparatus  is  disassembled  through  the  removal  of 
eight  bolts  and  reassembled  by  the  insertion  of  a  new  seal  and  bonded 
specimen.  While  some  of  the  seals  leaked  perhaps  due  to  plastic  flow  of 
the  copper  during  the  test,  it  is  felt  that  a  better  seal  design  might 
alleviate  this  minor  problem. 

before  a  new  test  is  begun,  the  assembly  is  bled  of  air  through  a 
small  valve;  it  is  then  placed  in  a  safety  cage.  While  this  test  does 
take  longer  to  complete  than  the  torsion  or  bending  experiments  of  Cl], 
it  can  be  done  in  about  five  minutes  by  an  experienced  operator. 


Pigure  4-2.  View  of  Hydraulic  Test  Assembly 


Valve  and  Coupling 


.Test  Fixture  for  Disk  Specimens 


Following  the;  Alumina  tests,  it  was  found  that  the  apparatus  could 
benefit  from  a  device  to  help  measure  incipient  specimen  cracking. 
Particularly  with  alumina,  the  pressure  at  failure  was  not  necessarily  the 
pressure  at  which  a  crack  first  appeared.  The  alumina  specimens  seemed 
to  arrest  cracks  growth  until  very  high  pressures  were  achieved.  A  conduct¬ 
ing  brittle  coating  over  the  outer  side  of  the  specimen  with  a  balance 
bridge  recorder  might  be  of  considerable  benefit  in  detecting  these  initial 
cracks.  Figure  4.8  depicts  such  devices. 


Specimen  gage  or 
brittle  coating 


Figure  4-8.  Specimens  Coating  Scheme  and  Bridge  Crack  Detector 


Chapter  5 


Results  of  the  Hydraulic  Tests 

Both  Alumina  (998)  and  Mullite  (MV33)  disks  made  from  the  same  batches 
as  the  tubular  specimens  of  [1]  were  tested  to  failure  through  hydraulic 
pressure  applied  to  one  surface.  The  experiment  was  designed  to  approximate 
simple  support  at  the  outside  where  the  soft  copper  seal  was  compressed 
between  the  specimen  and  pressure  retainer. 

Actual  test  results  performed  on  alumina  are  ranked  for  (as-cast) 
disks.  Table  5.1a,  and  for  ground-finish  disks.  Table  5.2a.  Four  of  these 
specimens  in  each  group,  Table  5.1b  and  5- 2b,  could  not  be  fractured  due 
to  leaky  seals. 

The  mean  strength  for  ground  specimens  over  'as-cast'  specimens 
(alumina)  showed  about  a  16.2%  increase  based  on  a  total  of  21  successful 
as-cast,  and  17  successful  ground-finish  disks. 

The  leaky  seals  encountered  in  both  sets  of  alumina  specimen  tests 
all  occured  at  nearly  no  pressure  reading.  Thus,  they  were  not  the  result 
of  normal  gasket  failure  in  the  copper  but  rather  a  result  of  potting 
material  failure  in  the  epoxy  which  was  employed  in  a  highly  unusual  fashion. 
Because  equal  numbers  of  as-cast  and  ground  finish  samples  could  not  be 
tested  due  to  leaks,  this  phenomenon  is  was  thought  not  to  be  related  to 
specimen  warpage. 

Of  further  concern  in  establishing  table  5- la  and  5-2a  is  the  fact 
that  the  actual  pressure  level  at  which  fracture  occured  in  the  alumina 
specimens  was  difficult  to  determine  with  precision.  Typically  specimens 
failed  without  an  apparent  noise  or  an  abrupt  pressure  drop.  Rather,  the 
alumina  samples  began  to  ooze  hydraulic  fluid  as  the  operator  attempted  to 
sustain  or  increase  pressure  levels.  The  figures  cited  in  Tables  5 -la  and 


Table  5 .la.  Alumina  Disk  (as  cast)  Pressure  Tests 
(Mean  Pressure  at  Failure  *  2860  p.s.i.) 


Specimen 

Number 

Applied  Pressure 
at  Failure  (p.s.i.) 

Remarks 

A24P 

1400 

All  tests  were  normal 

A25P 

1500 

A22P 

2000 

A10P 

2175 

A13P 

2225 

A14P 

2350 

A12P 

2400 

A1P 

2550 

A17P 

2750 

A16P 

2925 

A23P 

2970 

A26P 

3000 

A29P 

3050 

AGP 

3200 

A1P 

3200 

A11P 

3350 

A28P 

3400 

A15P 

2500 

A27P 

3670 

A18P 

3950 

A30P 

4500 

Table  5  .lb.  Alumina  Disk  (as  cast)  Pressure  Tests 
The  following  specimens  were  rejected. 


Specimen 

Number 

Maximum  Applied 

Pressure 

Remarks 

A3P 

Not  recorded 

Leakage  at  seal 

A4P 

Not  recorded 

Leakage  at  seal 

A5P 

Not  recorded 

Leakage  at  seal 

A6P 

Not  recorded 

Leakage  at  seal 

Table  5. 2a.  Alumina  Disk  (ground  finish)  Pressure  Tests 
(Mean  Pressure  at  Failure  ■  3365  p.s.i.) 


Specimen 

Nunfcer 

Applied  Pressure 
at  Failure 

Remarks 

A31P 

2360 

All  test  were  normal 

A41P 

2635 

A45P 

2750 

A35P 

2790 

A8P 

2800 

A7P 

2900 

A9P 

2955 

A42P 

3000 

A20P 

3150 

A3SP 

3150 

A33P 

3190 

A19P 

3650 

A39P 

3695 

A32P 

3985 

A40P 

4500 

A44P 

4600 

A21P 

5100 

Table  5. 2b.  Alumina  Disk  (ground  finish)  Pressure  Tests 
The  following  specimens  were  rejected. 


Specimen 

Number 

Maximum  Applied 
Pressure 

Remarks 

A34P 

Not  recorded 

Leakage  at  seal 

A36P 

Not  recorded 

Leakage  at  seal 

A37P 

Not  recorded 

Leakage  at  seal 

A43P 

3400  p.s.i. 

Produced  fracture  to  one  side 

5.2a  essentially  represent  a  pressure  at  which  significant  sustained 
leakage  occured  through  an  identifiable  crack.  It  generally  could  not 
be  exceeded  through  the  use  of  a  hand  operated  ram. 

The  test  results  for  Mullite  are  given  for  as-cast  samples  in 
Table  5.3a  and  for  ground  samples  in  Table  5.1*.  The  average  strength 
was  considerable  below  that  of  alumina  which  may  have  accounted  for 
fewer  failures  due  to  leakage.  The  as-cast  seal  failures  are  given  in 
Table  5. 3b.  There  were  no  failures  in  the  mullite  ground  finish  specimens. 
With  the  mullite  disk  test  there  is  a  10%  difference  between  average 
strengths  of  as-cast  tests  over  ground  test,  a  trend  that  is  opposite  that 
of  the  alumina  disks.  A  further  variance  of  the  mullite  results  is  that 
three  samples  in  the  Table  5.3b  failed  via  a  single  diagonal  crack.  These 
are  not  included  in  Table  5 -3a  due  to  the  concern  that  wrapage  of  the  as- 
cast  specimens  may  have  produced  these  anomalies. 

A  further  difference  in  the  testing  of  Alumina  and  Mullite  disks  is 
that  when  mullite  disks  fail,  there  is  an  instantaneous  pressure  drop  and 
an  unmistable  crack  in  the  specimen.  In  addition  to  the  tabulated  results 
of  the  pressure  test,  there  is  additional  information  to  be  gleaned  from 
an  observation  of  the  fractured  specimens.  In  the  photo  reproductions  that 
follow,  it  is  important  to  note  that  each  specimen  is  labeled  with  or  without 
a  trailing  "P"  in  the  specimen  number.  All  disk  specimens  in  Tables  5.1  - 
5.3  have  a  "P"  suffix  to  avoid  confusion  with  tube  specimens.  This  "P" 
indicates  the  specimen  side  exposed  to  hydraulic  fluid  pressure  and  the 
specimen  number  with  no  "P"  indicates  the  external  or  atmospheric  side. 

Thus,  for  hydraulic  specimen,  A13P,  a  photograph  of  its  sides  will  be 
labeled  A13  (outside)  and  A13P  (inside).  The  outside  surface  photos  can 
also  be  detected  by  the  adhesive  that  protrudes  at  the  edge. 


Table  5- 3a.  Mullite  Disk  (as  Cast)  Pressure  Tests 
(Mean  Pressure  at  Failure  =  1814  p.s.i,) 


Specimen  Number 

Applied  Presure 
at  Failure 

Remarks 

M120P 

1000 

M114P 

1100 

M119P 

1100 

M117P 

1200 

M121P 

1210 

M122P 

1290 

M123P 

1300 

M126P 

1300 

M109P 

1450 

M125P 

1450 

M101P 

1500 

K115P 

1500 

MllOP 

1750 

MI08P 

1850 

M103P 

1900 

M116P 

2100 

M113P 

2400 

M124P 

2450 

M99P 

2650 

M100P 

3700 

M102P 

3900 

Table  5. 3b.  Mullite  Disk  (as  Cast)  Pressure  Tests 
The  following  specimens  were  rejected. 


Specimen 

Number 

Maximum  Applied 
Pressure 

Remarks 

M104P 

M105P 

M106P 

M107P 

M111P 

M112P 

M118P 

950 

900 

800 

Not  Recorded 

Not  Recorded 

Not  Recorded 

Not  Recorded 

Single,  large,  diagonal  break 
Single,  large,  diagonal  break 
Single,  large,  diagonal  break 

Leak  at  seal 

Leak  at  seal 

Leak  at  seal 

Leak  at  seal 

52 


Table  5  .4.  Mullite  Disk  (Ground  Finish)  Pressure  Tests 
(Mean  Pressure  at  Failure  *  1638  p.i.l.) 
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A  series  of  two  alumina  photos  are  listed  for  both  specimens 


A11P  and  A13P.  Figure  5.1a  shows  the  outer  side  and  5  0b  shows  the 
inside.  The  inside  of  this  specimen  and  of  others  are  shown  in  greater 
magnification  to  provide  enhanced  evidence  of  cracks.  It  is  interest¬ 
ing  that  this  specimen  gives  no  indication  of  a  crack  on  its  compression 
or  inner  surface  although  there  is  a  radial  or  'crow  foot'  pattern  in 
the  outer  surface,  Figure  5.1a.  This  type  of  crack  was  typical  of  almost 
all  alumina  specimens  such  as  Figure  5. 2a  (outside)  and  5. 2b  (inside). 

While  the  outside  surface  exhibited  a  series  of  hairline  width  cracks,  no 
visual  crack  could  be  detected  on  the  inside.  Indeed  there  was  little 
change  in  pumping  pressure  as  detected  by  the  test  operator.  In  fact, 
with  the  alumina,  pressures  could  be  re-elevated  provided  the  pumping  rate 
exceeded  the  small  leakage  rate  through  the  specimen.  There  was  an  attempt 
to  detect  incipient  cracking  by  listening  while  observing  the  specimens  by 
closed  circuit  TV.  In  the  case  of  alumina,  this  revealed  nothing. 

Several  of  the  alumina  specimens  were  later  cleaned ;  they  appeared 
still  to  be  intact  except  for  the  hairline  cracks  made  visible  by  rubbing  the 
surface  with  graphite  mixture.  Again,  future  tests  could  benefit  from  an 
addition  of  a  conducting  brittle  strip  to  aid  in  crack  detection.  These 
specimens  appeared  to  have  cracks,  Figures  5.1  and  5.2,  that  propagated  to 
depth  of  about  50  mm  to  70  mm  under  the  outer  surface  where  they  appeared 
to  arrest.  The  exact  cause  of  the  hydraulic  leak  in  alumina  seems  to  have 
been  due  to  seepage  through  internal  cracks  that  became  sealed  after  the 
test  and  not  by  leakage  at  the  seal. 

The  corresponding  tests  on  mullite  showed  a  great  deal  of  difference 
in  both  the  strength  of  the  specimen  and  in  the  fact  that  again  fracture 
appeared  much  more  definite  and  catastrophic.  In  specimen  Ml 08,  repre¬ 
sentative  of  mullite  samples,  fracture  is  through  the  entire  thickness  of 
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the  specimen,  Figure  5.3a  (outside  surface)  and  Figure  5.3b  (inside 
surface).  While  the  fracture  is  through  the  entire  thickness,  the 
pattern  on  the  two  surfaces  is  greatly  different  indicating  that  the 
crack  has  not  propagated  uniformly  throughout  the  mullite.  There  is 
considerable  evidence  of  spalling  on  the  pressure  (compression)  side  of 
this  specimen  while  a  lesser  amount  of  material  has  been  lost  at  the 
outer  side. 

Specimen  M115  shown  in  Figures  5.^a  and  5.^b  provides  a  more  complex 
fracture  pattern  although  no  spalling  can  be  noted.  Again,  it  is  seen 
that  the  surface  crack  patterns  are  greatly  different  on  the  inside  and 
outside.  In  about  all  cases  the  radial  ray  or  'crow  foot'  fracture 
pattern  was  noted.  This  was  observed  in  both  alumina  and  mullite  failed 
samples  and  is  thought  to  be  the  result  of  a  similar  stress  field  with 
radial  symmetry  with  failure  occuring  along  rays  where  the  material  is 
statistically  weak. 


# 
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Figure  5.2a.  Hydraulic  Fracture  of  Specimens  A13 
Outer  Surface  (2.5X) 


Figure  5.2b.  Hydraulic  Fracture  of  Specimens  A13P 
Inner  Surface  (3X) 


Chapter  6 


Conclusion 

The  experimental  investigations  embodied  by  Navair  Reports 
N0001 9-78-C-0520  and  N000 19-79-0057^  form  a  set  of  experiments  that 
provide  Weioull  data  on  bending  and  torsion  for  tube  specimens  and  for 
a  uniformly  loaded  circular  plate.  General  conclusions  are: 

1 )  The  feasibility  of  torsion  and  bending  tests  on  one  specimen 
geometry  of  a  class  of  ceramic  material  has  been  demonstrated.  Evidence 
of  erratic  results  produced  by  stress  concentrations  near  points  of  load 
application  is  not  significant. 

2 )  The  use  of  inexpensive  tube  geometries  may  reduce  the  cost  of 
ceramic  specimen  tests  when  complemented  by  the  experiments,  theory 
and  programming  tecnniques  of  this  report. 

3)  The  use  of  hydraulic  pressure  in  providing  a  suitable  loading 
mechanisim  has  been  demonstrated.  Improvements  to  this  technique  nay 
however  be: 

a)  the  use  of  a  simple  retaining  bladder  to 
restrict  hydraulic  fluid  excursion  and  seepage 
into  potentially  porous  specimens. 

b)  the  inclusion  of  a  fracture  detecting  brittle 
conductive  coating  or  other  techniques  to  pro¬ 
vide  a  more  precise  measurement  of  incipient  fracture 
in  materials  that  have  some  fracutre  arresting 
capability  such  as  the  alumina  tested. 

c)  the  modification  of  current  copper  seals  and  specimen 
bonding  materials  or  techniques  to  reduce  reject  tests. 


4)  A  thorough  analysis  is  yet  to  be  made  of  the  disks  to  test 
whether  or  not  raw  Weibull  data  can  accurately  predict  failure  dis¬ 
tributions  in  materials  to  which  a  complex  stress  field  is  applied. 
Such  a  study  is  ongoing  and  will  be  the  subject  of  a  forthcoming 
submitted  publication. 
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A) 


Torsio  (AUTO. FOR) 


B) 


AGAR  (AGARWA. FOR) 


Appendix  II  Distribution  List  —  6? 


Q 


9 


9 


63 


program  torsto 

THIS  PPHCRAM  HIVES  WEIBULL  PARAMETERS  FOP  TORSION 
REAL  L 

DIMENSION  EPFf40).Wf 4n).AAf 3),  Pf40) 

DIMENSION  VAPf 3,3), PFSf 3,3) 

CnMMON/PAP/K.!,,C»T(40>#FJ,SlGU,EM,RO.PI,TPFf  40) 

WPITEf 5, 40 ) 

rOR“AT(Jx, 'ENTER  OAT A  FILE  NUMBER  AND  ITER  FRPQP  ROUND') 
*BNO  OF  DATA  POINTS 

J»TOROllF  AT  rupture 

PEAO(*>,*)  NO, FPP 
PEAOfNU,#)  K 
DO  400  Tc) ,K 
RFADfND,*)  TfT),WfT) 

„  T f  T )sl 2 ,*T f  T ) 

THIS  IS  Tn  CONVERT  FROM  FT. LB,  TO  IN.I.R, 

THE  FILE'S  UNITS  APF  ET.IP. 


This  IS  Tn  CONVERT  FROM  FT.LB,  TO  IN.I.R, 

V  THE  FILE'S  UNITS  APF  FT.LP. 

*0V  CONTINUF 

,  XsFLOAT f  K4 1  ) 

t  NOW  RANK  THE  DISTRIBUTION 
,  DO  1  !■!  , K 

*  EPF  f  I ' *FLD AT  f  I )  /X 
readInd.#)  L,RO,RT 

_  ,  E.T*3. 1  41  A/2 RO*«4-RT#*4) 

c  L«HAGF  LENGTH  OF  SPFCTHrN 
t  KOHOllTrR  RADIUS  OF  SPECIMEN 
C  KlsINNFB  RADIUS  OF  SPECIMEN 
C  LjaPOLAP  MOMENT  OF  SPECIMEN 
-  WRTTEC5.59)  L.PO.PI.E.t 

?9  EORMATC/.7X,  •!.«'. IX, FB. 4,  •  P0»  '  ,  1 X  .  Ffe  .  4f  •  PT* • , 1 X , FR . 4 . 

1  •  EI«'f 1X,E10.S) 

...  GO  TO  II 

}ll  WPTTEf 5. 112)  C,Fm,SIGH#PES(1,3) 

*14  FORMAT ( IX, 'C*FM, SIGU,RFS* ' , 4E1 S ,6) 

*  1  WPTTEfS.42) 

♦2  FORMAT! IX, 'INPUT  FM,STGU,  AND  ICODE  FOR  C  ITERATION') 

PEADf 5,#)  EM, SIGU, ICODE 
IE  f ICODE, EO.l)  GO  TO  90 
8TGMAXRTCK)*PO/Ei) 

DENI «(STGMAX»5IGU)*#EM 
CNC«ALQGf2,)/DENl 
VAPf 1 , 3)«ALOGf 2,*FLOAT(K) 1/DENl 
_  CN«0 

C  NOW  BEGIN  C  ITERATION  TO  MINIMISE  RES 
WRTTEf5,69)  VARf),3),ENC 

09  roRMATdx, 'input  initial  c  and  dft.ta  c,ol"  c,d«»,?eio,4) 

C  PFADf 5**)  VAP( 1,1) ,ENC 

*0  CO«CN 

f  2  ICTrO 

*4  VAPf l,l)BVAPf 1 .2) 

VAP(l,2)BVAPfl,|) 

RFS  f 1 , i )aPE3( 1 #2) 

RESf l,2)aPE6f)*3) 

VARf  1,1)bVAPM,1)4ENC 
ICT«ICTt1 

IE  f ICf.GT.lOOO)  GO  TO  98 
ravisM 


.rfe.4,'  PI* ' , 1 X , FR , 4 , 


CsVAR f i , J I 
CALL  A&AR 

now  compute  The  r 

RES f 1 , ; )>0 , 
00  S  !■) ,N 
WNafEPf f I)* 
RES! 1 , 1)«RE 


RESf 
00  S 
WNaf 
RES! 


RESIDUALS 


)»TPFf I)) 

RESf l,l)4WfI)#WNBWN 


W 
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c 

/8 

C  If 

C  lr 
t  Rn 


c  ir 


WRTTE(5,7t)  rvtPM  ,.71  >,  J1>1,3).  (RES  ft  ,  J?),J2«i ,3) 
EORMATflX, 'VAR.RES  J  TD  3  IRE!  'AE10.4) 

T.ESS  THIN  3  VALUES  JUST  CONTINUE, 


^  u  u  "  a  i  i.u-i  ■  *iwb| 

IE  f ICT.LT. 4 )  GO  TO  14 

RESM.2)  IR  T.ESS  THAN  RE.S(1,1)  ANn  RESfl.3)  THEN  WE  HAVE  A 
CALCULATE  WHERE  MTN  OCCURS,  DROP  ENC  BY  10  AND  CONTINUE. 

TE  C(RES(1,2).LT.RESf  1  ,!)). AND. {RFSC  1,21. T.T, REST  1,3M) 

1  GO  TO  13 

REST  1,3 ) .LE.RFSf 1 »  2  )  CONTINUE  ELSE  REVERSE  ENC, 

IF  TRESfl  ,3)  .LE.RESM  ,2))  GO  TO  14 
ENC«-ENC/2. 

GO  TO  12 

CALL  MXNIf 1 ,VAR,RES,CN) 

WRiTE(5,RR)  CN.RESfl.33 

FORMAT  (IX,*  CN  ANP  RESM.3)  FQP  THIS  MTN  ARE*  »  ,  2E1  5.0) 

IF  (ARSf (CO-CN)/CN).LT.ERR)  GO  Tn  111 
ENC*EwC/10. 

VAR(I,3)«CN 


MIN. 


GO  TO  ID 
DO  81  lal.K 

NRTTE(5f78l  fT,FPEfn,TPF(I),WfTn 
FORMAT (IX,  •  I, EPE.TPF.WW ',12.3(1  X.E1  3.M) 
continue 

WRITE(5.310)  „  „ 

FORMAT! IX, *  NOV  INPUT  SlGl .DELSIG1 .NO  OF  OIVISIONS') 
REAP ( 5 , * )  ST1.DST1.N00 
DO  31?  I«1.N0P*1 

F*0 

IF  (STl.LE.STGU)  GO  TO  313 

Z«C*(ST1-SIGU)#*EM 

F«1 ,-FXPf-Z) 

WRTTE(S| 314)  ST1.F. 

FORMAT! 1 X, *  STRESS* ' , El  0  #  5, *  PR0R«',F7,5) 

ST1«ST1*DST1 

CONTINUE 

WRITE!®, 88)  ICT 

format!?*, 'IMAX  EXCEEDED,  IMAX«'.T6) 

GO  TO  99 

STOP 

END 


SUBROUTINE  WJNI(I, VAR. RKB.XXX) 

DIMENSION  VARf3«3).RES(3. 1)  . 

PI aRESl T,1)/((VAR(I,1 )»VAR(I,2) )#(VAR(I, 1 )»VAP (I , 3) ) ) 
P2«BES(I,2)/((VAR(I,2)*VAR(I,1))#(VAR(I,2)»VAR(I,3))) 
P3«RESCI.3)/((VAR(I.31-VAR(I,1 ) )* ( VAR ( I , 3 ) -VAR ( I , 2 ) ) ) 
AAANP1TP24P3 

RBRwPl Rf VAR(7 , 2)4VAR(T,3))*P2*(VAR(I»1 ) ♦VAR(I , 3) ) 
♦P3#(VARf f ,1)+VAP( J,2)) 

XXX«.5*BBR/AAA 

IE  !VARfI,n.GT,VAR(I,3))  GO  TO  2 
IE  ((XXX.GT.VAR(I,3)).0R,(XXX.LT.VAR(I,1))) 
XXX«(VAR(I,1 )+VAR(I ,3))/2, 

RETURN 

TE  (fXXX.0T,VARfT*l)),0P»(XXX,LT.VARfI»3))) 
XXX«(VaR(I, 13+VARf T,3>)/2, 

RETURN 

END 


XXX) 


VAR (I , 3) ) ) 
VAR (I , 3) ) ) 
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SUBROUTINE  AGAR 
REAL  L 

r  _  eOMMON/PAP/K,|.,C,T(4<n.FJ,SIGU,FM,RO.PI,TPFf 40) 
c  Compute  determined  ppobablilties  at  each  topoue, 

J;  &TOPE  TN  TPF(1  TO  K>. 

L  t IPST  COMPUTE  NECESSARY  CONSTANTS, 

B»3.1 416*2. «L»C/fFM*l.) 

DO  1  I  *  t , K 

rprm»o. 

,,  TTL*TfI)/FJ 

t  NON  COMPUTE  R  AT  WHICH  THE  THRESHOLD  STRESS  OCCURS, 
RLOWSRI 
PTH«STGf»/TII, 

IF  (RTH.GF.PO)  GO  TO  1 
^  IF  fRTH.GT.Rn  PLOWbRTH 

C  Now  COMPUTE  BN. 

C1«P0*TIL 

C7*(C1«SIGU^**(EM4lt) 

C3«G2*(C1-SIGU) 

DI«PLnW*TIL 
D?»(ni«SIGin*«(FM+J  .) 

D3*02*(01-S!GU) 

BV*P*f (BO*C2»C3/fTTL*fEM+2.))3- 
1  (PLOW*U2»D3/fTIL«(EM^?.))))/TlL 

.  TPF(T)M1 ,-EXPf-PN) 

1  CONTINUE 

RETURN 
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